Determination of significant immunological timescales from mRNA-LNP-based vaccines in humans.

A compartment model for an in-host liquid nanoparticle delivered mRNA vaccine is presented. Through non-dimensionalisation, five timescales are identified that dictate the lifetime of the vaccine in-host: decay of interferon gamma, antibody priming, autocatalytic growth, antibody peak and decay, and interleukin cessation. Through asymptotic analysis we are able to obtain semi-analytical solutions in each of the time regimes which allows us to predict maximal concentrations and better understand parameter dependence in the model. We compare our model to 22 data sets for the BNT162b2 and mRNA-1273 mRNA vaccines demonstrating good agreement. Using our analysis, we estimate the values for each of the five timescales in each data set and predict maximal concentrations of plasma B-cells, antibody, and interleukin. Through our comparison, we do not observe any discernible differences between vaccine candidates and sex. However, we do identify an age dependence, specifically that vaccine activation takes longer and that peak antibody occurs sooner in patients aged 55 and greater.

Multiple cohort study of hospitalized SARS-CoV-2 in-host infection dynamics: Parameter estimates, identifiability, sensitivity and the eclipse phase profile.

Within-host SARS-CoV-2 modelling studies have been published throughout the COVID-19 pandemic. These studies contain highly variable numbers of individuals and capture varying timescales of pathogen dynamics; some studies capture the time of disease onset, the peak viral load and subsequent heterogeneity in clearance dynamics across individuals, while others capture late-time post-peak dynamics. In this study, we curate multiple previously published SARS-CoV-2 viral load data sets, fit these data with a consistent modelling approach, and estimate the variability of in-host parameters including the basic reproduction number, R0, as well as the best-fit eclipse phase profile. We find that fitted dynamics can be highly variable across data sets, and highly variable within data sets, particularly when key components of the dynamic trajectories (e.g. peak viral load) are not represented in the data. Further, we investigated the role of the eclipse phase time distribution in fitting SARS-CoV-2 viral load data. By varying the shape parameter of an Erlang distribution, we demonstrate that models with either no eclipse phase, or with an exponentially-distributed eclipse phase, offer significantly worse fits to these data, whereas models with less dispersion around the mean eclipse time (shape parameter two or more) offered the best fits to the available data across all data sets used in this work. This manuscript was submitted as part of a theme issue on "Modelling COVID-19 and Preparedness for Future Pandemics".

School and community reopening during the COVID-19 pandemic: a mathematical modelling study.

Operating schools safely during the COVID-19 pandemic requires a balance between health risks and the need for in-person learning. Using demographic and epidemiological data between 31 July and 23 November 2020 from Toronto, Canada, we developed a compartmental transmission model with age, household and setting structure to study the impact of schools reopening in September 2020. The model simulates transmission in the home, community and schools, accounting for differences in infectiousness between adults and children, and accounting for work-from-home and virtual learning. While we found a slight increase in infections among adults (2.2%) and children (4.5%) within the first eight weeks of school reopening, transmission in schools was not the key driver of the virus resurgence in autumn 2020. Rather, it was community spread that determined the outbreak trajectory, primarily due to increases in contact rates among adults in the community after school reopening. Analyses of cross-infection among households, communities and schools revealed that home transmission is crucial for epidemic progression and safely operating schools, while the degree of in-person attendance has a larger impact than other control measures in schools. This study suggests that safe school reopening requires the strict maintenance of public health measures in the community.

Long-term durability of immune responses to the BNT162b2 and mRNA-1273 vaccines based on dosage, age and sex.

The lipid nanoparticle (LNP)-formulated mRNA vaccines BNT162b2 and mRNA-1273 are a widely adopted multi vaccination public health strategy to manage the COVID-19 pandemic. Clinical trial data has described the immunogenicity of the vaccine, albeit within a limited study time frame. Here, we use a within-host mathematical model for LNP-formulated mRNA vaccines, informed by available clinical trial data from 2020 to September 2021, to project a longer term understanding of immunity as a function of vaccine type, dosage amount, age, and sex. We estimate that two standard doses of either mRNA-1273 or BNT162b2, with dosage times separated by the company-mandated intervals, results in individuals losing more than 99% humoral immunity relative to peak immunity by 8 months following the second dose. We predict that within an 8 month period following dose two (corresponding to the original CDC time-frame for administration of a third dose), there exists a period of time longer than 1 month where an individual has lost more than 99% humoral immunity relative to peak immunity, regardless of which vaccine was administered. We further find that age has a strong influence in maintaining humoral immunity; by 8 months following dose two we predict that individuals aged 18-55 have a four-fold humoral advantage compared to aged 56-70 and 70+ individuals. We find that sex has little effect on the immune response and long-term IgG counts. Finally, we find that humoral immunity generated from two low doses of mRNA-1273 decays at a substantially slower rate relative to peak immunity gained compared to two standard doses of either mRNA-1273 or BNT162b2. Our predictions highlight the importance of the recommended third booster dose in order to maintain elevated levels of antibodies.

Mathematical modelling of vaccination rollout and NPIs lifting on COVID-19 transmission with VOC: a case study in Toronto, Canada.

BACKGROUND: Since December 2020, public health agencies have implemented a variety of vaccination strategies to curb the spread of SARS-CoV-2, along with pre-existing Nonpharmaceutical Interventions (NPIs). Initial strategies focused on vaccinating the elderly to prevent hospitalizations and deaths, but with vaccines becoming available to the broader population, it became important to determine the optimal strategy to enable the safe lifting of NPIs while avoiding virus resurgence. METHODS: We extended the classic deterministic SIR compartmental disease-transmission model to simulate the lifting of NPIs under different vaccine rollout scenarios. Using case and vaccination data from Toronto, Canada between December 28, 2020, and May 19, 2021, we estimated transmission throughout past stages of NPI escalation/relaxation to compare the impact of lifting NPIs on different dates on cases, hospitalizations, and deaths, given varying degrees of vaccine coverages by 20-year age groups, accounting for waning immunity. RESULTS: We found that, once coverage among the elderly is high enough (80% with at least one dose), the main age groups to target are 20-39 and 40-59 years, wherein first-dose coverage of at least 70% by mid-June 2021 is needed to minimize the possibility of resurgence if NPIs are to be lifted in the summer. While a resurgence was observed for every scenario of NPI lifting, we also found that under an optimistic vaccination coverage (70% coverage by mid-June, along with postponing reopening from August 2021 to September 2021) can reduce case counts and severe outcomes by roughly 57% by December 31, 2021. CONCLUSIONS: Our results suggest that focusing the vaccination strategy on the working-age population can curb the spread of SARS-CoV-2. However, even with high vaccination coverage in adults, increasing contacts and easing protective personal behaviours is not advisable since a resurgence is expected to occur, especially with an earlier reopening.

Efficacy of a "stay-at-home" policy on SARS-CoV-2 transmission in Toronto, Canada: a mathematical modelling study.

BACKGROUND: Globally, nonpharmaceutical interventions for COVID-19, including stay-at-home policies, limitations on gatherings and closure of public spaces, are being lifted. We explored the effect of lifting a stay-at-home policy on virus resurgence under different conditions. METHODS: Using confirmed case data from Toronto, Canada, between Feb. 24 and June 24, 2020, we ran a compartmental model with household structure to simulate the impact of the stay-at-home policy considering different levels of compliance. We estimated threshold values for the maximum number of contacts, probability of transmission and testing rates required for the safe reopening of the community. RESULTS: After the implementation of the stay-at-home policy, the contact rate outside the household fell by 39% (from 11.58 daily contacts to 7.11). The effective reproductive number decreased from 3.56 (95% confidence interval [CI] 3.02-4.14) on Mar. 12 to 0.84 (95% CI 0.79-0.89) on May 6. Strong adherence to stay-at-home policies appeared to prevent SARS-CoV-2 resurgence, but extending the duration of stay-at-home policies beyond 2 months had little added effect on cumulative cases (25 958 for 65 days of a stay-at-home policy and 23 461 for 95 days, by July 2, 2020) and deaths (1404 for 65 days and 1353 for 95 days). To avoid a resurgence, the average number of contacts per person per day should be kept below 9, with strict nonpharmaceutical interventions in place. INTERPRETATION: Our study demonstrates that the stay-at-home policy implemented in Toronto in March 2020 had a substantial impact on mitigating the spread of SARS-CoV-2. In the context of the early pandemic, before the emergence of variants of concern, reopening schools and workplaces was possible only with other nonpharmaceutical interventions in place.

Analysis of Host Immunological Response of Adenovirus-Based COVID-19 Vaccines.

During the SARS-CoV-2 global pandemic, several vaccines, including mRNA and adenovirus vector approaches, have received emergency or full approval. However, supply chain logistics have hampered global vaccine delivery, which is impacting mass vaccination strategies. Recent studies have identified different strategies for vaccine dose administration so that supply constraints issues are diminished. These include increasing the time between consecutive doses in a two-dose vaccine regimen and reducing the dosage of the second dose. We consider both of these strategies in a mathematical modeling study of a non-replicating viral vector adenovirus vaccine in this work. We investigate the impact of different prime-boost strategies by quantifying their effects on immunological outcomes based on simple system of ordinary differential equations. The boost dose is administered either at a standard dose (SD) of 1000 or at a low dose (LD) of 500 or 250 vaccine particles. Results show dose-dependent immune response activity. Our model predictions show that by stretching the prime-boost interval to 18 or 20, in an SD/SD or SD/LD regimen, the minimum promoted antibody (Nab) response will be comparable with the neutralizing antibody level measured in COVID-19 recovered patients. Results also show that the minimum stimulated antibody in SD/SD regimen is identical with the high level observed in clinical trial data. We conclude that an SD/LD regimen may provide protective capacity, which will allow for conservation of vaccine doses.